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Abstract 

We analyze the structure of the perturbation expansion of the general multi- 
channel Kondo model with channel anisotropic exchange couplings and in the 
presence of an external magnetic field, generalizing to this case the Anderson- 
Yuval technique. For two channels, we are able to map the Kondo model 
onto a generalized resonant level model. Limiting cases in which the equiv- 
alent resonant level model is solvable are identified. The solution correctly 
captures the properties of the two channel Kondo model, and also allows an 
analytic description of the cross-over from the non Fermi liquid to the Fermi 
liquid behavior caused by the channel anisotropy. 



I. INTRODUCTION 

The single channel Kondo model has a long history as the simplest model believed to 
contain the relevant physics of magnetic impurities embedded in metals. A lot of efforts 
have been devoted to study this model and presently one can safely claim that it has been 
completely understood from the theoretical point of view. Apart from the original per- 
turbative scaling approach, which already gave the correct qualitative descriptionSi, there 
is also an exact solution of this model obtained by Bethe-ansatz techniqueil. The physics 
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underlying the single channel Kondo model is the formation of a non degenerate singlet 
state at low temperature. The impurity spin is screened by the conduction electrons, hence 
the magnetic susceptibility, obeying Curie- Weiss law at high temperatures, undergoes the 
Kondo cross-over saturating to a constant upon lowering the temperature. 

In an attempt to describe realistic magnetic impurities (which have orbital structure), 
different generalizations of the simple Kondo model have been proposed^. Among these 
generalizations, the simplest one is probably the model describing an impurity spin S coupled 
to channels of conduction electrons (commonly referred to as the multichannel Kondo 
model). Surprisingly, it has been realised that, for > 2S, this model exhibits a behavior 
qualitatively different from the single channel onei, which is due to a non-trivial ground 
state with a residual degeneracy. This gives rise to divergent low temperature susceptibility 
Ximp and specific heat coefficient 7 = Cv,imp/T (so-called non Fermi liquid behavior). It 
should be said that the experiments on dilute magnetic alloys (the systems, for which the 
Kondo model was originally proposed) do not give any clear evidence for such a behavior but 
rather suggest that the ground state is always a singlet. A possible explanation would be 
that the channel symmetry (which is necessary for the non Fermi liquid behavior to occur) 
is broken (since typically no exact symmetry guarantees channel equivalence). If the energy 
scale of this symmetry breaking term is small, a slow cross-over from the non Fermi liquid 
to the Fermi liquid behavior is expected. However, it might be difficult to experimentally 
distinguish it from the usual Kondo cross-over. Consequently, since early 80's, the interest 
to the issue was decaying, even though, in the meanwhile, some exact solutions for the case 
of equivalent channels became available. 

In recent years, however, interest to the subject grew again, as other, more promising, 
realizations of the multichannel Kondo model have been proposed; for instance, two level 
systems in metalsi, heavy fermion compounds^, high Tc superconductorsi. While for the two 
level systems there are convincing experimental evidences^ of the non Fermi liquid behavior 
predicted by the theory, the other proposed realizations are till now controversial. 

The general multichannel Kondo hamiltonian is given by 



N N 



(3) 



Hk = EE^o(V^..,0 + E {j^aS^'ym + ^ [5V-(0) + 5-a+(0)]| , (1) 

a=l cr a=l ^ ^ ^ 

where 

i/o(^/',V^^) = Eek4^k (2) 

k 

is the kinetic energy of the conduction electrons ipaa, a- = 1, ■■■,N being the channel index, 
and a =t, i being the spin index (we assume a spin 1/2 impurity). Notice however that 
in some realizations of this model, the channel index is the physical spin while the spin 
index labels an orbital quantum number (see also Section IV). The electron spin densities 
in Eq. (|^) are defined by 

<W = I [^l^(x)^„^(x) -^^^(x)V'„^(> 

In the single channel case (A^ = 1) the exchange couplings J flow to infinity under scaling, 
if positive. The spin anisotropy 7^ J± disappears at the fixed point. This is interpreted 
as the formation of a singlet at the impurity site. On the contrary for the > 1 channel 
symmetric case, the infinite coupling fixed point is unstable as well as the weak coupling 
one, hence a stable intermediate coupling fixed point. 

Until recently, the information on the behavior of the model around this fixed point 
could only be extracted from the Bethe-ansatz solution0 and conformal field theorylHI. In 
1992, Emery and Kivelsonlll provided a simple solution for the symmetric two channel case 
at a particular value of the longitudinal exchange coupling. The solution was obtained 
by means of abelian bosonization technique generalizing a procedure previously applied to 
the single channel Kondo model by Schlottmannlli. The original Schlottmann's approach 
was in turn inspired by Toulouse's mapping of the single channel Kondo model onto a 
resonant level model, which he achieved by the analysis of the partition function^, borrowing 
the perturbative treatment of Anderson and YuvalS. The two methods, bosonization and 
Anderson- Yuval method, are equivalent, in the sense that they give qualitatively similar 
results (it has been checked in the single channel case). Still, the latter is formally more 



rigorous (e.g. it does not fully relies on band linearization) and straightforward (it simply 
amounts to compare the perturbation expansions of two different models). 

In this paper we generalize the Anderson- Yuval method to the multichannel Kondo model 
with channel anisotropic exchange couplings. For the two channel case, we are able to show 
that the perturbation expansion of the Kondo model is equivalent to that of a generalized 
resonant level model. In the channel isotropic case, the model is of the type found by Emery 
and Kivelson via bosonization technique. The novel feature is the channel anisotropy which 
gives rise to interesting cross-over phenomena. We demonstrate that even in this case a 
mapping of the Kondo model (for particular values of the longitudinal exchange couplings) 
onto a solvable resonant level model does exist. For > 2 we are unable to find any kind 
of resonant level model which would reproduce the perturbation expansion. 

II. GENERALIZATION OF THE ANDERSON YUVAL APPROACH 

In this Section, we analyze the structure of the perturbation expansion of the hamiltonian 
(|^) in the transverse exchange couplings. 

A. Single channel model 

Consider first the single channel problem. We allow for an exchange anisotropy [J^ = 
Jy = J± is different from J^). We want to calculate the impurity partition function in 
time space, using a perturbation expansion in powers of J_l. A term of order 2n involves 
2n alternate impurity spin flips. Let tj be the times of up flips, t'- that of down flips (n 
of each). The philosophy is to calculate that particular term exactly, for a given history 
and to show that it is identical to the corresponding expansion for another problem 
(a resonant level), with appropriately chosen parameters. The two problems are mapped on 
each other term by term: they are equivalent. Note that we thus bypass summation of the 
perturbation series in Jj_. For a specific value of (the "Toulouse limit") the equivalent 
problem happens to be trivially solvable: if we can scale through that value we have an 
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explicit description of the crossover to low temperature - the " (100 — e)%" exact solution of 
Anderson. The error stems from the fact that universal scaling is not just a change of J^. 
That error is supposed not to change the qualitative behavior - and anyway it is implicit in 
the equivalent bosonization technique (less powerful since it relies on a Born approximation 
for phase shifts). 

Assume first that Jz = 0. Each vertex flips a conduction electron spin. A ti vertex 
creates a | electron and an | hole, a t- vertex does the reverse. In a typical diagram, the 
electron propagators go from any ti to any t- for | spins, and from any t[ to any ti for f 
spins. Since = these propagators are free electron local propagators 

G„(i) 



t - i^o ^signt' 

where z/q is the density of states for one spin at the Fermi level, and is a high energy 
cut-off of the order of the conduction bandwidth. Let D^iti^t^) be the contribution of 
the n lines with spin a that join ti and t[ vertices. D^j has a pole whenever ti = t'j (one 
propagator is divergent). Moreover crossing symmetry implies a zero whenever ti = tj 
(or t- = t'j): exchanging the extremities of two propagators changes the sign. Hence the 
Cauchy determinant found by Anderson and Yuval (here and in what follows we omit time 
independent pre-factors containing the density of states uq and to the short-time cut-off, 
restoring them in the final expressions) 

D = '-^ (A) 

hj 

Expression (H) is homogeneous with degree (— n), as expected for D„. The proof is completed 
by looking at the asymptotic behavior. The corresponding contribution to the impurity 
propagator is Uq = D^D^. 

We now restore J^. The potential felt by a spin a electron changes at each impurity 



flip - hence an edge singularity that will modify the long time behavior of U , see Ref. |T6 



In Fig.l we have drawn the time dependent potential felt e.g. by the up-spin conduction 



electrons at the impurity site (the down-spin one is the opposite) The phase shift is 5+ (resp. 
5_) when the electron and impurity have parallel (resp. antiparallel) spins. What matters 
is the discontinuity of phase shifts when a flip occurs, 5 = 5+ — If we assume electron 
hole symmetry, then 

5+ = = tan"^(7rz/oJ^/4) (5) 
The effect of the flip is two-fold: 

{{) The open lines that contribute to D^j can scatter any number of times on the impu- 
rity. That generates once again the Cauchy determinant (via a Muskhelishvili type 
of analysi^l). For each spin D^^ is replaced by /)i-25/7r^ thus open line contribution 
acquires an extra factor Ul = D~^^^'^. 

[ii) In addition the underlying Fermi sea reacts to the flip via closed loops that exponentiate 
(for a given history the potential is structureless). The resulting contribution is Uc = 



[Hi) Altogether U = UqUlUc, and, after inserting back the prefactors, we obtain 

en 
SO 



where rj is an exponent that depends on J; 

^ = 2 f 1 - - 



6V 



TT 



(Note that rj = if 6 = n, which corresponds to the strong coupling limit 5+ = — 5_ = 
n/2). 

Let us now consider a resonant level model for spinless electrons, characterized by the 
hamiltonian 



V 



H = Ho{^,^^) + X ^t(o)rf + c^t^(o) +- ^^(0)^(0) -^(O)^t(O) idU--) (7) 



2 
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where d is an impurity orbital at the Fermi energy located at the origin. \1/ is a free Fermi 
field which kinetic energy is the same as in Eq.(|l|). The interaction potential V produces a 
phase shift discontinuity 

,' = 2ta„-'(^^) (8) 

between the empty and full d-states. We expand in powers of A which plays the role of J±. 
The ti and t[ vertices correspond to d'^ and d operators. The structure of the expansion is 
the same as for the Kondo case, except that there is no spin degeneracy. We have one set 
of open lines originating from A vertices, which can scatter off the flipping S'. We also have 
one set of closed loops - hence altogether 



2n 



D 



U'=(X^^o^o) I^J (9) 



with 

v'=(i-- 



5'V 



TT 



Thus the Kondo problem with coupling 6 is mapped term by term onto the resonant 
level with coupling 6' if the two propagators Eqs .(^ and (|^) are identical. This implies 
1] = 7]' (that can always be achieved by appropriately choosing V) and 



The Toulouse limit corresponds to 5' = 0, i.e. a phase shift 5 = 7r(l — 1/^/2). The resonant 
level hamiltonian can then be trivially diagonalized, yielding a low temperature Fermi liquid 
behavior. The resonant level hamiltonian of the form has been previously derived by 
Wiegmann and Finkelstienlli. 



B. Multichannel model 



We now turn to the N channel case. In a first stage we assume fiavor degeneracy: what 
does remain of the previous analysis? In order to answer that question we proceed in reverse. 
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(i) The alternation of up and down spins is unchanged. The flipping potential due to 5 
is the same as before (see Fig.l), whatever the flavor involved at each vertex. Edge 
singularities are consequently unaffected. The scattering contribution to open lines 
is again Ul = D"^^^'^ (flavor is fixed by extremities). Each closed loop can have an 
arbitrary fiavor and therefore Uc = D^'^^'^ /'^'^ , 

(ii) Paradoxically the difficulties come from the part Uq (in the absence of J^). Uq still 
has poles whenever = t'j, but one looses crossing symmetry. If the ends of two lines 
are interchanged one usually changes the number of closed loops C - hence a change 
in the degeneracy N'-' . As a result Uq cannot be expressed simply in terms of D. 

In order to proceed, we must assign to each vertex its fiavour index a. The and t'j 
then break into subclasses, tia and t'j^ (a = 1, ...,N). For a given diagram the number 
of vertices in different subclasses needs not be the same, but spin and flavor conservation 
implies an equal number of t^a and t'j^ within a given subclass. Since flavor is conserved 
along an open line, Uq is a product of independent factors Uoa- For each factor crossing 
symmetry holds and f/oa is the square of a Cauchy determinant Da as it would be for a 
single channel. Da is still given by (§), the products running only over a-type times, t^a and 
t'ja- In the end we flnd (again omitting prefactors proportional to J±, uq and ^o) 

U = [D^...Dr,]' D-'-^^'^'l^T = l^lz^D^N (10) 



P^ = 1-4- + 2n(-] . (11) 



where 

2 

Note that D is not the product of individual Da', we have instead D = Di...DnF, in which 
we have set 



n n (^»» ~ ^jb) [t'ia - t'jf, 

ij a<b 
ij a<b 



The factor F couples the channels. 



The Emery-Kivelson solution to the two channel case, N = 2, is based on a mapping of 
the Kondo problem onto the following spinless resonant level hamiltonian 



;M/t(0)M/,(0) - vl/,(0)M/t(0)] {d^d - 1) (12) 



V 
+ 2 



where d is again a fictitious spinless Fermi operator. Notice that we have introduced two 
Fermi fields \1/ and coupled to the impurity in a different way (the reason why we have 
not used the same field will become clear later). In order to establish the equivalence we 
first consider the case V = We divide the d and d'^ flips into two subclasses, depending on 
whether a fermion \1/ is emitted or absorbed. Times tn and tj2 correspond to flips with 
a fermion emitted or absorbed respectively, t'^^ and t-2 are their hermitian conjugates. A 
spinless fermion propagator can go as usual from tn to t'^^ (t[2 to tj2) or from tn to tj2 (t^2 
to tj-^). The latter possibility is the new feature. The corresponding impurity propagator 
U' has poles whenever a propagator has zero time range, i.e. when tia = t'j^, tn = tj2 or 
t-2 = t'j^. Due to crossing symmetry it has zeroes when tia = ija, i'ia = ^'ja^ = ^'j2 
^'ii ~ Once again one thereby builds a Cauchy determinant which happens to be 

U' = ^ (13) 

Expression (O) has the right poles and zeroes. It moreover has the right overall power of t 



and asymptotic behavior: it is the correct answer. Comparing (p!3| ) with the definition of F 
we see that 



D 

We now restore the flipping potential V. Since it involves a different Fermi field, it gives 
rise only to a closed loop contribution. Altogether we have 



U' = D\^) (14) 

with the same 5' as in Eq. (^. Comparing (0) with (|TU|) we see that the two problems are 
mapped onto each other if [see Eq. (0)] 
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which is always possible since both right and left sides are positive. Notice that for a given 
6 the interaction potential V in Eq. (|12]) is, according to (||) and (|15D , given by 



V=—ts^n(-] = Atanf--^) (16) 

The problem is directly solvable if V = 0, i.e. when 6 = tt/2. In the electron hole 
symmetric case that implies 5+ = — (5_ = 7r/4 - a typical intermediate coupling as expected 
for the two channel overscreened Kondo impurity. Indeed from the expression of j32 we see 
that the model is symmetric under 6 ^ n — 6. This extends the result of Ref. ^ that the 
two channel Kondo model behaves similarly around Jz = (i.e. 6 = 0) and Jz = oo (i.e. 
(5 = tt). By symmetry the fixed point should be exactly at (5 = 7r/2, that is at the solvable 
line V = 0. 

The argument can be extended to a flavor dependent exchange J. Due to anisotropy we 
must treat separately the channel dependence of Jj_ and Jz- Different Jj_i and J±2 do not 
affect the structure of the perturbation expansion. As we have shown above the mapping 
works as follows: 

^l^(0,tii)^it(0,ta)5+(tii) ^ S{ta)^{0,ta) 
vl/l^(0,4)vI/,^(0,4)S-(4) ^ vl;t(0,4)rf(4) 

thus we need only modify accordingly the flipping matrix elements of the equivalent model, 
which becomes 



J±i 
2 



J 



uo^o [rf^^(O) + ^\0)d\ + -j^^jMo [d^^\0) + ^{0)d\ (17) 

Notice that, if J±i = J ±2, ( P^ reduces to (|^ with A = J±Wi^o^o- 

A difference between Jzi and Jz2 gives rise to different phase shifts 61 and 62. Let us first 
consider the scattering correction to the "1" open line, Uli- The Muskhelishvili propagator 
for a channel 1 spin up electron is 
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t-t 



it - Ua) it' - t'J 
it - 4) it' - Ua) 



Its contribution to Uli is obtained by putting t equal to any tn, t' to any t-^, hence a factor 
square it in order to account for spin and we multiply by the corresponding 
term for channel 2. The closed line contribution is straightforward since flavor is conserved 
along a loop, 

Uc = 



Altogether the impurity propagator is 

JJ ^ _Q2-45i/7r^2-452/7r^-2(5i+52)/7r£)2(<52+52)/7r2 

(Remember that F = D / D1D2 ) If 5i = 82 = 5, we recover the previous result ([T0|). 
In general let us write 

5i = 5 + e , 82 = 5 — e 



(18) 



Then (pTSf) reduces to 



U 



D \Di 

The additional factors with respect to ( p^ can be reproduced with an extra potential 



(19) 



W 



*^(0)^(0) -^(0)^^(0)1 (dU 



With the choice 



W 



2 2 [81-62 

tane = tan 

TTZ/Q vrz/Q V 2 ^ 



(20) 



the closed loop contribution generates the factor _D(2e/7r) rjj^^ j^^^ factor of ( [T9| ) comes 
from the fact that the phase shift discontinuity is +e on the "1" vertices, —e on the "2" 
vertices (see Fig. 2). In that way one can map any version of the two channel Kondo impurity 
onto an extended Emery-Kivelson hamiltonian. 
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It is interesting to examine to what extent such an analysis could be pursued if > 2. 
We return to the flavor symmetric case, for which (|^ holds. If we manage to have Pn = 0, 
then 

_ iD,...DN? Di-Dn 
^ - D ~ F ^^-^^ 

Hence two questions: (z) Can we achieve /3 = 0? (m) If we can, is there a solvable model 
that gives the same Ul It is clear that no real phase shift 5 will achieve /? = if > 2. 
That may be a definitive objection since poor man's scaling scans the real 5 axis. Let us 
ignore it, hoping that some analytic continuation argument might help. Then in order to 
reproduce (^Tj) we must introduce a coupling 

^+(^1 + ... + ^^) +H.C. 

in which the operators are such that the corresponding propagators are 

{g{t) ~ 1/t is the free electron propagator). Then U will have poles whenever tia = t'^^ on 
the one hand, tia = tjb, = t'^^ (a ^ h) on the other. It will have zeroes if tia = ^ja, t'^^ = i'ja 
or tia = tjf,. That just generates the combination (0). It remains to be seen what kind of 
algebra could produce (p^): we do not know of any. 

III. MAGNETIC FIELD EFFECTS 

Let us consider the effects of an uniform magnetic field B = (0, 0, B) in the framework of 
the Anderson- Yuval approach. The magnetic field appears in the hamiltonian with a term 

N 



Hb = —I^bB 



giS"" + gcj rfx^o-a^(x 

a=l 



(23) 



where the electron spin density is defined in Eq. (^, gi and gc are the Lande factors of the 
impurity and the conduction electrons respectively, and /is is the Bohr magneton. 
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As before we will treat the transverse exchange perturbatively. This implies that the 
reference states |t) and |J.), which are used for the perturbation expansion, are the eigenstates 
of the hamiltonian with fixed impurity spin direction in the presence of the magnetic field 

N J N 

Hyi = E ^a.) ± E -f^aM T ^5 - g^^^sB ^ / rfxa„,(x). (24) 

We have to understand how the magnetic field modifies the perturbation expansion in 
J^. B gives rise to two effects. 

{i) It shifts the chemical potential for up and down spin electrons (in opposite directions). 
This causes a small change in the spin up and down phase shifts if the band has a 
finite curvature at the Fermi energy. This effect is negligible at low temperature. 

[a) It causes a difference = E-^ — in the ground state energies of (|24D for the two 
impurity spin directions, which appears in the Muskhelishvili propagators. 

By standard phase shift arguments, based on Friedel's sum rule for the displaced charge, 
we find 

a=l-^^F~ 2 

where ep is the Fermi energy. For small magnetic field (|25| ) reduces to 



1 ^ r^F+^B 

AE = -fisg^B + - E / .... dM^) (25) 



AE = -^iB9^B + J2 ^a. (26) 

^ a=l 

The above energy difference enters in the impurity propagator ([T0|) via the following phase 
factor 



exp 



n 

-iAEY.{U-t[) 



(27) 



i=l 

The conduction electron part of AE actually represents the leading term of closed loops 
diagrams, that one which grows linearly with (tj — t[) instead of logarithmically. 

Which term has to be added to the resonant level model in order to reproduce (P^? It 
is easy to realize that the corresponding term is simply 
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AES, I — > AE iSd--] . (28) 



2 

Notice that at the Emery-Kivelson hne for the two channel case 



TT 



a=l,2 

SO that if Qi = Qc then AE = (at first order in B). Consequently the impurity mag- 
netic susceptibility vanishes, in agreement with conformal field theoryS and bosonization 
approacheJii. This in turns means that at the Emery-Kivelson line the reference states are 
such as to perfectly screen the impurity spin. When the departure away from the Emery- 
Kivelson line is treated as a perturbationEl, both the specific heat and susceptibility acquire 
logarithmic sigularities, leading to the universal Wilson ration Ryy = 8/3. 

As to the N > 2 channel symmetric case, conformal field theory0 and abelian be 



josoniza- 



tion approaches (which till now exist only for = 4, see Ref. ^Oj) again predict the impurity 
susceptibility to vanish at the fixed point. From Eq. (p6D , we see that AE = for 6 = n/N, 
and therefore the impurity susceptibility is rigorously zero. If this is the true property of 
the fixed point, as it follows from the analysis of Ref. ITl|, then 6 = tt/N is the fixed point. 



IV. SOLUTION OF THE TWO CHANNEL ANISOTROPIC MODEL 

In this section we discuss the two channel Kondo model in more detail, focusing on the 
effects of channel anisotropy. The starting hamiltonian is 

= E E ^o(^.., O + E {j^aS'<iO) + [5V-(0) + S-a^iO)] I . (29) 

a=l cr a=l ^ ^ J 

If the exchange couplings are channel symmetric Jzi = Jz2 and J±i = J±2, it is knownll 
that the hamiltonian ( |29| ) flows towards a non trivial fixed point. At this fixed point the 
model exhibits non Fermi liquid behavior, namely the impurity susceptibility Ximp and the 
specific heat over temperature Cv,imp/T diverge at low temperatures as ln(l/T), and the zero 
temperature entropy is finite and equal to ln(2)/2, as if half of the impurity spin degrees of 
freedom were decoupled from the conduction electrons. Physically this occurs because two 
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(and more) channels tend to overscreen the impurity spin, so that the complete screening 
characteristic of the single channel Kondo model can not take place, thus leaving a ground 
state degeneracy. In the case of a finite channel anisotropy, the system will always choose 
the channel with the strongest exchange to screen the impurity spin, and the usual Fermi 
liquid behavior of the single channel model will finally take place at zero temperature. The 
corresponding RG flow diagram is sketched in Fig.3. 

As we have shown in the previous sections, this hamiltonian (|^) can be mapped onto 
the following resonant level hamiltonian 



+ 



J 



±2 



+ 



W 



(0)^(0) -^(0)^"f(0)l (dU 



+ 



V 



vI/t(0)vI/,(0)-vl>^(0)vI/t(0)l 



1 



AE ( dU 



(30) 



2/ V 2. 

where the interaction potentials are related to the longitudinal exchange couplings via [see 

Eqs.o-(|ra)-(mi 



V 



Jzl ~ J; 



z2 



2 1 + 7rVoV,iJ,2/16' 
8 l-nXJziJz2/16 



and (assuming equal impurity and conduction electron Lande factors gi = Qc 



(31) 
(32) 



AE 



-tan 



TT 



V 2 



j B = -gfiBHV)B. 



In the case of symmetric exchange couplings, the resonant level model ( pOD reduces to 
the hamiltonian (|1^) originally considered by Emery and Kivelson. Then the combination 
d'^ + d is decoupled from the conduction electrons (hence the ground state degeneracy and 
the non Fermi liquid behavior). A finite channel anisotropic transverse exchange couples 
this combination to conduction electrons and moves the system away from the non Fermi 
liquid fixed point towards the Fermi liquid single channel fixed point0. The smaller is the 
anisotropy, the lower is the cross-over temperature. In what follows we analytically study 
this cross-over in the solvable limit V = W = 0. 
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Since the total number of fermions is not conserved by the hamiltonian, there are anoma- 
lous Green functions. In the Nambu representation 



D 



the impurity Green function 



G,{t) = -i{T (D{t)D\0))) 



is a 2 X 2 matrix. Its Fourier transform can easily be evaluated. For u much smaller than 
the bandwidth, we find: 



To " Tx 



+ 



To + 



2 uj + iFsignu; 2 u + 27signa; 
where the resonance widths are defined by 



(33) 



vr 



r = ^^^0^0 (^±1 + J±2f , 

7 = J^O^O {J±i - J±2f , 

fi being the Pauli matrices, and Tq the unit matrix. The impurity spectral function is 

F 1,. . . 7 



and it is therefore equally shared by two lorentzians with different widths F and 7. In the 
channel isotropic case 7 — 0, one of the two lorentzians tends to S{uj), representing the 
impurity degree of freedom which is decoupled from the conduction band in this particular 
limit0. 

The impurity contribution to the free energy can be calculated in a standard way by 
integration over the coupling constant. The result is 



F(T)=Fo(T) + 1^/(0; 



tan { — \ + tan ^ (— 



(34) 



where F^iT^ is the free energy in absence of coupling between the impurity and conduction 
electrons, jiijS) is the Fermi distribution function and the integral should be limited to the 
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conduction bandwidth. The entropy can be calculated by S{T) = —dF(T)/dT. By defining 
the function 



S{z) 



2nz 



1 



2'KZ 



1 



Inr ( - + — 

,2 2-nz 



Invr 



(35) 



where V'(^) is the psi-function and T{z) is the gamma-function, the entropy turns out to be 

T < 7 



S{T) = ln(2) + s(^l'^+S 



( ttT 1 V 



In V2 , 7 < T < r 

r + 7 



(36) 



In 2 



27rT 



T > r 



the last equality being valid for 7 <^ F. S(T) is shown in Fig. 4. We see that S{0) = 0, as 
expected since no degeneracy is left, but there is a region of temperatures (the wider the 
smaller 7 is) where the entropy is close to that of the symmetric two channel model. 

Another quantity of physical interest is the longitudinal impurity susceptibility. As we 
know from the above analysis, exactly on the Emery-Kivelson line, Ximp — has 
to consider deviations from this line (i.e. 7^ 0) in order to account for a finite impurity 
susceptibilityEl. The resulting susceptibility is 



imp 



Jo 



[9f^BAiV)Y 



1 



1 F 

^ ' 2 ^ 27rT 



^ \2 2ttT 



7r(F-7) 

In the case 7 <^ F the susceptibility shows the same kind of cross-over behavior as the 
entropy: 

1 



X 



imp 



l9teA(V)|'' 



7r(F-7) V7 

7r(F-7) VT 
1 

I 4T ' 



T < 7 
7 < T < F 
T > F 



(37) 



As expected the magnetic susceptibility saturates at low temperature, although at interme- 
diate temperatures it shows the logarithmic behavior of the two-channel Kondo model. 

It follows from ( PB] ) and ( P7| ) that the Wilson ration Rw is not universal: it depends on 
the amount of anisotropy, 7/F. Such a conclusion is obvious in the limit of small anisotropy, 
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when the energy scales are well separated. Then the residual entropy In \f2 must be quenched 
in a temperature range ~ 7, implying Cy^imp ~ T j^, while the susceptibility Ximp just rounds 
off logarithmic singularity, Ximp ~ lii(r/7): the Wilson ratio is very small. Such a lack of 
universality is also apparent in the phenomenological, Fermi liquid description of the low 
temperature limit, T <^ 7. Then the impurity is quenched into a singlet, and the residual 
conduction electron phase shift in the channel (m, cr) may be expanded as 

where Srim'a' is the change in the occupation measured from the ground state. Universality 
implies that Sm(j{e) is invariant 

{{) If the chemical potential of the other channel is changed (there is no channel flip) 

(m) If e and the chemical potential are changed by the same amount (the Kondo singularity 
is attached to the Fermi level) 

Hence in our two channel case 



^2<7 = S20 + 0!2 



e — 



+ 6iaa'5n2a' 
+ 92aa'5nia' 



(the cross terms 6*1, 6*2 are equal in the electron- hole symmetric case 610 = 620 = tt/2.) Us is 
the one channel density of s-states at the Fermi level. It is then straightforward to extend 
the analysis of Ref. ^: the resulting impurity corrections are 

Cy^imp _ ai + a2 Ximp _ 2(ai + 0^2) ^ O1 + O2 

Cy TTUs ' X TTi/g 27r 

Due to the channel interaction 9 the Wilson ratio Rw = CyXimp/ xCv,imp departs from the 
single channel value 2. Put another way, one has a line of fixed points rather than a unique 
one. If " 1" is the dominant screening channel, Ji goes to infinity while J2 may evolve towards 
any arbitrary value: once the spin S is screened, J2 no longer scales. This arbitrariness is 
reflected in the Wilson ratio. 



To our knowledge, the two channel Kondo model is most convincingly realized by two 
level systems in metal alloys!. This has recently been experimentally confirmed thanks 
to the development of point contact spectroscopy^. In these systems, the role of the spin 
is played by some orbital degree of freedom, while the physical spin plays the role of the 
channel index. Thus the model is by construction channel isotropic. However, an external 
magnetic field breaks the channel symmetry and generates an effective channel anisotropy 
proportional to the curvature of the conduction electron band times the magnetic field B. 
In this case, the coupling to the magnetic field is described by the following term in the 
hamiltonian: 

where a is now the pseudo-spin index, and the channel indices 1 and 2 correspond to the 
physical spin up and down respectively. The magnetic field shifts the Fermi level of channel 
1 with respect to that of channel 2. Electron-hole symmetry is thereby broken within each 
channel. It follows that the magnetic field induces a phase shift anisotropy, 6i — 62 oc B, 
reflected into a finite W oc B. Such a correction comes both from the correction to the 
Fermi level density of states 

e=tF 

and from the Zeeman shift of band edges. The change of the density of states also modifies 
the pseudospin flip amplitudes in the equivalent resonant level model: 



2 



2 



The hybridization anisotropy is equivalent to a finite 7 



7 = ^A9^^B)%'-^JlB\ (38) 

Thus, B causes the cross-over to a Fermi liquid behavior at low temperature as observed in 
Ref. ^ As to the physical magnetic susceptibility, it is related to the first derivative of the 
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free energy ( p4D with respect to 7. The low temperature (low magnetic field) behavior of 
the susceptibility is given by: 

Ximp — 77 <iO'J± in 



mm 



T' fi2 



16z/o 

(One can show, that finite W oc B does not contribute to the log-divergent part of the 
susceptibility.) 

Very recently, the channel anisotropic (but spin isotropic) Kondo model has been solved 
using Bethe ansatz methods by N. Andrei and A. JereJi. Their conclusions are qualitatively 
similar ours. 

We are thankful to N. Andrei for helpful discussions. 
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FIGURES 

FIG. 1. Time-dependent potential seen by a spin up conduction electron. The impurity spin 

flips from | to j at times ti and vice versa at times t'j. 

FIG. 2. Time-dcpcndcnt local potential felt by an electron in the effective resonant level model. 
While the same annihilation operator ^ is involved in the two types of flips, the potentials are 
opposite. The flips at times t'^^ and tj2 are the hermitian conjugates. 

FIG. 3. Qualitative (Ji, J2) RG flow diagram for the anisotropic two channel Kondo model. 

FIG. 4. Entropy -S'(r) for various values of the anisotropy = 7/r: from the top 
A = 0,0.1, 0.5,1. 
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